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CORRECTION TO BLACK-SCHOLES FORMULA DUE TO 
FRACTIONAL STOCHASTIC VOLATILITY 


JOSSELIN GARNIER*AND KNUT S0LNAt 

Abstract. Empirical studies show that the volatility may exhibit correlations that decay as a 
fractional power of the time offset. The paper presents a rigorous analysis for the case when the 
stationary stochastic volatility model is constructed in terms of a fractional Ornstein Uhlenbeck 
process to have such correlations. It is shown how the associated implied volatility has a term 
structure that is a function of maturity to a fractional power. 
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1. Introduction. Our aim in this paper is to provide a framework for analysis of 
stochastic volatility problems in the context when the volatility process possesses cor¬ 
relations that decays like a power law. We will both consider the case of “long-range” 
processes where the consecutive increments of the process are positively correlated, 
corresponding to the so called Hurst coefficient H > 1/2, as well as the case with 
“short-range” processes with consecutive increments being negatively correlated with 
H < 1/2. Replacing the constant volatility of the Black-Scholes model with a random 
process gives price modifications in financial contracts. It is important to understand 
the qualitative behavior of such price modifications for a (class of) stochastic volatility 
models since this can be used for calibration purposes. Typically the price modifica¬ 
tions are parameterized by the implied volatility relative to the Black-Scholes model 
[27, 42]. For illustration we consider here European option pricing and then the im¬ 
plied volatility depends on the moneyness, the ratio between the strike price and the 
current price, moreover, the time to maturity. The term and moneyness structure 
of the implied volatility can be calibrated with respect to liquid contracts and then 
used for pricing of related but less liquid contracts. Much of the work on stochastic 
volatility models have focussed on situations when the volatility process is a Markov 
process, commonly some sort of a jump diffusion process. However, a number of 
empirical studies suggest that the volatility process possesses long- and short-range 
dependence, that is the correlation function of the volatility process has decay that is 
a fractional power of the time offset. This is the class of volatility models we consider 
here. We find that such correlations indeed reflect themself in an implied volatility 
fractional term structure. An important aspect of the modeling is also the presence of 
correlation between the volatility shocks and the shocks (driving Brownian motion) of 
the underlying, this “leverage effect” influences the implied volatility in an important 
way and we shall include it below. The leverage effect is well motivated from the 
modeling viewpoint and important to incorporate to fit observed implied volatilities, 
albeit a challenging quantity to estimate [2]. Evidence of leverage and persistence 
or long-range dependencies have been found by considering high-frequency data and 
incorporated in discrete time series models [8, 20, 43]. 

Here we model in terms of a continuous time stochastic volatility model that is 
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a smooth function of a Gaussian process. We use a martingale method approach 
which exploits the fact that the discounted price process is a (local) martingale. We 
model the fractional stochastic volatility (fSV) as a smooth function of a fractional 
Ornstein-Uhlenbeck (fOU) process. We moreover assume that the fSV model has 
relatively small fluctuations, of magnitude (5^1 and we derive the associated leading 
order expression for implied volatility with respect to this parameter via an asymptotic 
analysis. This gives a parsimonious parameterization of the implied volatility which 
may be exploited for robust calibration. The fOU process is a classic model for 
a stationary process with a fractional correlation structure. This process can be 
expressed in terms of an integral of a fractional Brownian motion (fBm) process. 
The distribution of a fBm process is characterized in terms of the Hurst exponent 
H G (0,1). The fBm process is locally Holder continuous of exponent a for all a < If 
and this property is inherited by the fOU process. The fBm process, IV/^, is also 
self-similar in that 

R} *=■ {a^W/^,tG K} for all a > 0. (1.1) 

The self-similarity property is inherited approximately by the fOU process on scales 
smaller than the mean reversion time of the fOU process that we will denote by 1/a 
below. In this sense we may refer to the fOU process as a multiscale process on 
relatively short scales. 

The case with H G (0,1/2) gives a fOU process that is a so-called “short-range” 
dependent process that is rough on short scales and whose correlations for small time 
offsets decay faster than the linear decay associated with a Markov process. In fact the 
decay is as the offset to the fractional power 2H. In this regime consecutive increments 
of the fBm process are negatively correlated giving a rough process also referred to 
as an anti-persistent process. The enhanced negative correlation with smaller Hurst 
exponent gives a relatively rougher process. 

The case with H G (1/2,1) gives a fOU process that is a so-called “long-range” 
dependent process whose correlations for large time offsets decays as the offset to the 
fractional power 2{H — 1). It follows that the correlation function of the fOU process 
is not integrable. This regime corresponds to a persistent process where consecu¬ 
tive increments of the fBm are positively correlated. The relatively stronger positive 
correlation for the consecutive increments of the associated fBm process with increas¬ 
ing H values gives a relatively smoother process whose correlations decay relatively 
slowly. For more details regarding the fBm and fOU processes we refer repectively to 
[7, 17, 18, 37] and [10, 35]. 

In order to simplify notation and interpretation of the results we present them in 
the context of the fractional OU process. However, as we show in Appendix B, the 
results readily generalize to the case with general Gaussian processes with short- and 
long-range dependence. 

A large number of recent papers have considered modeling of volatility in terms of 
processes with short- and long-range dependence. In [13] the authors consider a long 
memory extension of the Heston [34] option pricing model, a fractionally integrated 
square root process, a generalization of the early work in [14]. They make use of 
the analytical tractability of this model, in fact a fractionally integrated version of 
a Markovian affine diffusion, with affine diffusions considered in [19]. The emphasis 
is on the long-range dependent case {H > 1/2) and long time to maturity. The au¬ 
thors focus on the conditional expectation of the integrated square volatility and show 
the fractional decay of this, moreover, they discuss estimation schemes for model pa¬ 
rameters based on discrete observations. In the Markovian case the mean integrated 
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square volatility would exponentially fast approach its mean value and flatten the 
implied volatility term structure. They remark that long-range dependence provides 
an explanation for observations of non-flat term structure in the regime of large matu¬ 
rities since the long-range dependence may make the implied volatility smile strongly 
maturity-dependent in this regime, while also producing consistent smiles for short 
maturities. The model presented in [13] was recently revisited in [31] where short and 
long maturity asymptotics are analyzed using large deviations principles. 

The concept of RFSV, Rough Fractional Stochastic Volatility, is put forward in 
[5, 30]. Here a model with log-volatility modeled by a fBm is motivated by analysis 
of market data, which they state provide strong support for a value for the Hurst 
exponent H around 0.1. As explained above small values for H correspond to very 
rough processes. It is remarked that such a process can be motivated by modeling of 
order flow using Hawkes processes. The authors discuss issues related to change from 
physical to pricing measure and use simulated prices to fit well the implied volatility 
surface in the case of SPX with few parameters. They argue that the fractional model 
generates strong skews or “smiles” in the implied volatility even for very short time 
to maturity so that this modeling provides an alternative to using jumps to model 
such an effect. The form of the implied volatility surface and the structure of the 
returns have been used to argue that the asset price should be a jump process [1, 9]. 
Indeed models with jumps may be used as an alternative approach to capture smile 
dynamics to the fractional approach considered here and recent contributions consider 
models driven by Levy processes both for volatility models [21, 40] and directly for 
price models [4]. 

A variant of the model in [30] is considered in [39] where the log-volatility is 
modeled as a fractional noise, with fractional noise being the increment process of a 
fBm for a certain increment length. The authors discuss the well-posedness of this 
model from the financial perspective and in doing so make use of a truncated version 
of the integral representation of the fBm. In [38] this model is supported by data 
analysis and motivated by an agent-based interpretation. 

In [11, 12] the authors consider the situation when the volatility is modeled as a 
function of a fOU process whose shocks are independent from those of the underlying. 
Their focus is on a tree-based method for computing prices, estimation schemes for 
model parameters, and a particle filtering technique for the unobserved volatility given 
discrete observations. They consider some real data examples and find estimated 
values for the Hurst exponent which is larger than 1/2, in particular in a period after 
a market crash. In [32] small maturity asymptotic results are presented for this model. 

Among the many papers considering short maturity asymptotics, in the early 
paper [3] Aids et al. use Malliavin calculus to get expressions for the implied volatility 
in the regime of small maturity. They find that the implied volatility diverges in the 
short-range dependent case and flattens in the long-range dependent case in the limit 
of small maturity. These results are consistent with what we present below. The 
modeling in [3] differs from the modeling below in that the authors consider volatility 
fluctuations at the order one level while below the fluctuations are relatively small, 
however, we consider any time to maturity. 

Fukasawa [28] discusses the case with small volatility fluctuations and short- and 
long-range dependence impact on the implied volatility as an application of the general 
theory he sets forth. He uses a non-stationary “planar” fBm as the volatility factor 
so that the leading implied volatility surface is identified conditioned on the present 
value of the implied volatility factor only, while below with a stationary model the 
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surface depends on the path of the volatility factor until the present, reflecting the 
non-Markovian nature of fBm. In [29] Fukasawa discusses the case of short-range 
dependent processes and short time to maturity and a framework for expansion of the 
implied volatility surface. He uses a representation of fBm due to Muralev [41]. He 
also considers local stochastic volatility models and find that these are not consistent 
with power laws in this regime. 

As a further generalization relative to a fractional Brownian motion based model 
the case of multi fractional Brownian motion based models is considered in [16] . This 
allows for a non-stationary local regularity or a time dependent Hurst exponent and 
then the implied volatility depends on weighted averages of the local Hurst exponent. 

In [23] Forde and Zhang use large deviation principles to compute the short ma¬ 
turity asymptotic form of the implied volatility. They consider the correlated case 
with leverage and obtain results that are consistent with those in [3]. They consider 
a stochastic volatility model based on fBm and also more general ones where the 
volatility process is driven by fBms and which is analyzed using rough path theory. 
They also consider large time asymptotics for some fractional processes. 

Indeed, a number of recent papers have considered small maturity asymptotics 
for implied volatility in the context of mixing, short- or long-range processes. Many of 
these use large deviation principles or heat kernel expansions [6, 23, 33], while another 
approach is to consider the regime around the money [3, 29, 40]. Recent works deal 
also with the regime of large strikes and derive bounds on the implied volatility [36] . 
Here we take another approach by considering a perturbation situation so that the 
implied volatility can be expanded around an effective volatility [27], also for large 
times to maturity. We model the volatility as a stationary process, a continuous 
time stationary short- or long-range dependent stochastic volatility process, with a 
view toward constructing a time consistent scheme. We use an approach based on 
the martingale method which is adapted to the fact that the volatility process is 
not a Markov process. We explicitly take into account the effects of correlation in 
between volatility shocks and shocks in the underlying, the leverage effect, and its 
form in short- and long-range dependent cases. We obtain expressions for the implied 
volatility for all times to maturity and also for log-moneyness of order one. Explicitly, 
we model the volatility as 

at=a + F{6Zi^), ( 1 . 2 ) 

for the fOU process that we discuss in more detail in Section 2.2. The function 
F is assumed to be one-to-one, smooth, bounded from below by a constant larger 
than —CT, with bounded derivatives, and such that E(0) = 0 and F'{0) = 1. It 
follows that the volatility process inherits (qualitatively) the correlation properties of 
the fBm process. Indeed, we have 

at=a + F{5Z^) =a + SZ^^ + 5^h\Z^) (1.3) 

where h^{y) = {F{Sy) — Sy)/d^ can be bounded uniformly in S by: 

\h^{y)\ < \\F''\\ooy^. (1.4) 

Note that throughout the paper we will be working with non-dimensionalized quan¬ 
tities. Specifically, if t' represents dimensionalized time say in units of “trading year” 
and T' is a typical time horizon being for instance a typical maturity time in years 
then t is the non-dimensionalized time: 

t=^- ( 1 - 5 ) 
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The main result is then the associated form for the implied volatility, see Equa¬ 
tions (5.1), (5.3) and (5.4) below, we summarize the result next. The implied volatility 
is here the volatility value that needs to be used in the constant volatility Black-Scholes 
European option pricing formula in order to replicate the asymptotic fSV option price, 
it is, up to terms of order 5^: 


for 


It 


= E 


r 1 
-T-t 



+ A{T -1) 


1-b 


\og{K/Xt)i 

{T-t)/fV 


( 1 . 6 ) 


A{t) 


6paT^^2 j-^ 

2r(i7-b |) r 




(1.7) 


where 1/a is the mean reversion time of the fOU process and r = 2a characteristic 
diffusion time for the underlying. Furthermore, Xt is the underlying price process with 
evolution as in (3.1) and Tt its associated filtration. Moreover, p is the correlation 
in between the Brownian motions driving respectively the volatility process and the 
underlying price process, K is the strike price so that KjX is the moneyness, and 
finally t = T — t is time to maturity. The first term in the implied volatility is the 
expected effective volatility over the remaining time period of the option conditioned 
on the knowledge at time t, note that this term is random. The second term is a 
leverage term which is present in the case that the underlying and the volatility have 
correlated evolutions so that p is non-zero. Note that p is commonly assumed to be 
negative. The log-moneyness term becomes relatively more important as the time to 
maturity becomes small relative to the characteristic diffusion time. 

In the short and long time to maturity regimes we then have for the leverage 
term: 
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We moreover have for the predicted effective volatility term: 
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for ar <C 1, 
for ar ^ 1, 
( 1 . 8 ) 

(1.9) 


( 1 . 10 ) 


It is important to note that we only assume r = T — t>0 so that in fact the implied 
volatility for small times to maturity may be very large for short-range dependent 
processes. This reflects the fact that for short-range dependent processes the volatility 
path is rough and may have a significant impact beyond the current predicted effective 
volatility level. However, when used in the standard Black-Scholes pricing formula 
the implied volatility indeed gives a pricing correction that is 0(5) for any r > 0. We 
also note that in the long maturity regime the implied volatility level may diverge for 
long-range dependent processes reflecting the fact that long-range dependence gives 
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strong temporal coherence and therefore relatively large corrections to the predicted 
current effective volatility. 

Note next that the calibration of the leverage component of the implied volatility 
in the general case in (1.6) involves estimation of the group market parameters: 

CT, H, (Sp), a, (1.11) 


from observed implied volatility data. In order to fully identify the model at the 
current time t we need moreover to estimate the current predicted effective volatility 
over a time to maturity horizon, that is, <Jt,t+T for 0 < r < T^ax ~ i- 

It is important to note that in our framework the market parameters are from 
the theoretical point of view independent of the current time t. Thus, in order to 
calibrate the model with data over a current time epoch ti < t < t2 one may use all 
the implied volatility recording in a joint fitting procedure. 

We remark that our results would be modified under the presence of general inter¬ 
est rates and market price of risk factors that we do not consider here. We also remark 
that identifying a “smile” shape, that is a more general function in log-moneyness, 
would require a higher-order approximation of implied volatility [26]. Finally, observe 
that the case H = 1/2 corresponds neither to a short-range dependent process nor a 
long-range dependent process, but the standard case of an Ornstein-Uhlenbeck pro¬ 
cess and a stochastic volatility that is a Markovian process with correlations decaying 
exponentially fast [27]. 

The framework we have presented is general and can be used for processes for 
which we can identify the key quantities of interest below. We discuss one important 
special case corresponding to a slow fOU process. In this case we model the volatility 
in terms of the “slow” fOU process 



whose natural time scale is 1/5 and whose variance is order one and given by CTou 
defined by (2.5) below, independently of 5. Then the volatility is 

at = (1.13) 


where F’ is a smooth, positive-valued function, bounded away from zero, with bounded 
derivatives. We introduce the two parameters 

ao = F{Zl'^), po = F'{Z^,-^), (1.14) 


that is, the local level and rate of change of the volatility. In this case the implied 
volatility is given by: 
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for To = 2,/a/^. Thus, the slow fractional volatility factor yields an implied volatility 
that corresponds to the one of the fractional model in (1.2) in the regime of small 
maturity, as given in (1.8). In the special case that FI = 1/2 the volatility process 
becomes a standard Ornstein-Uhlenbeck process and is in the class of slow processes 
considered in [27] and indeed the implied volatility in (1.15) can then be show to be 
exactly of the form discussed for the slow correction in [27] (Chapter 5). 
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The outline of the paper is as follows. First in Section 2 we introduce the details 
of the ingredients of the fSV model. In Section 3 we derive the main result of the 
paper, the leading order expression for the price in the situation with a fSV. The 
derivation is based on a contract with a smooth payoff function while the European 
payoff function has a kink singularity and we generalize the result to this situation 
in Section 4. Then in Section 5 we derive the expression for the implied volatility 
and how the fractional character of the volatility affects this. We connect to the slow 
time volatility model in Section 6 and present some concluding remarks in Section 7. 
In Appendix A we characterize some quantities of interest and associated technical 
lemmas that are being used in the price derivation in Section 3. 

2. The fractional stochastic volatility model. We describe in more detail 
the fBm and fOU processes that are used in the fSV construction (1.2). 

2.1. Fractional Brownian motion and its moving-average stochastic in¬ 
tegral representation. A fractional Brownian motion (fBm) is a zero-mean Gaus¬ 
sian process (IF/^)tgR with the covariance 

^ (If" + -\t- sp"), (2.1) 

where an is a positive constant. 

We use the following moving-average stochastic integral representation of the fBm 

[37]: 
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where {Wt)tGR is a standard Brownian motion over R. In this model (lF")tgR is a 
zero-mean Gaussian process with the covariance (2.1) where 
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2.2. The fractional Ornstein-Uhlenbeck process. We then introduce the 
fractional Ornstein-Uhlenbeck process (fOU) as 
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( 2 . 6 ) 












Note that it is not a martingale, neither a Markov process. 

Substituting (2.2) into the second representation in Eq. (2.4) gives in view of 
stochastic Fubini the moving-average integral representation of the fOU: 


where 
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f Kit — s)dWs 

J —CO 
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( 2 . 8 ) 


The properties of the kernel K are the following ones: 

- K is nonnegative-valued, K G L^(0,oo) for any H G (0,1) with /“/C^(it)du = o-q, 
and K G L^(0,oo) for any H G (0,1/2). 

- For small times at <C 1: 
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r(i7+l)a^-2 


■((at)^ 2+o((at)^ 2 )). 
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- For large times at ^ 1: 

K{t) = 


r(i7-i) 
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( 2 . 10 ) 


For H G (0,1/2) the fOU process possesses short-range correlation properties: 

E[Zf = cr^^^l - y(2H + o((as)2^)), as < 1. 

For H G (1/2,1) it possesses long-range correlation properties: 

+ o{{asf^-'^)y as > 1. 


( 2 . 11 ) 


( 2 . 12 ) 


The expansion (2.12) is valid for any H G (0,1/2) U (1/2,1) and for H G (1/2,1) it 
shows that the correlation function is not integrable at infinity. This is in contrast to 
the case of short-range dependent processes and also to Markov processes for which 
the correlation function is integrable. 

3. The option price. The price of the risky asset follows the stochastic differ¬ 
ential equation: 


dXt = atXtdW;, (3.1) 

where the stochastic volatility is 

at = d + F{6Z^), (3.2) 

has been introduced in the previous section and is adapted to the Brownian 
motion Wt, and IF/ is a Brownian motion that is correlated to the stochastic volatility 
through 


w; = pWt + 


( 3 . 3 ) 










where the Brownian motion Bf is independent of Wt- We remark that the main aspect 
of the model whose consequences we want to analyze here are the short- respectively 
long-range properties of the correlation function in Eqs. (2.11) and (2.12) under the 
presence of leverage as in Eq. (3.3). We will find that this has a dramatic effect on 
the asymptotic prices and the associated implied volatility. 

The function F is assumed to be one-to-one, smooth, bounded from below by a 
constant larger than —a, bounded above, with bounded derivatives, and such that 
F{0) = 0 and F'{0) = 1. Note that with this normalization for the function F 
it will not appear explicitly in the price approximation in Proposition 3.1 below as 
further properties are not important in that context. Moreover, then the filtration Ft 
generated by {Bt, Wt) is also the one generated by Xt. Indeed, it is equivalent to the 
one generated by (Wj*, Wj), or {Wt^Z^). Since F is one-to-one, it is equivalent to 
the one generated by (Wt*,crt). Since a + F is positive-valued, it is equivalent to the 
one generated by (W^*, cr^), or Xt- 

We aim at computing the option price defined as the martingale 

Mt=E[h{XT)\Ft], (3.4) 

where ft, is a smooth function with bounded derivatives apart from a finite set of 
points where it may have a jump discontinuity in its derivative. Note that the proof 
in this section will be given for the case with ft smooth and bounded. However, as we 
only need to control the function q[^\x) defined below rather than ft the argument 
can be extended from the smooth case to the situation with jump discontinuity in 
the derivative which is relevant in the case with a call payoff. We carry out this 
generalization explicitly in Section 4. 

The idea of the proof that we present below is to construct an approximation for 
Mt which has the correct terminal condition and which up to small (order 6'^) terms 
is a martingale. It then follows that we have a price approximation to 0{6‘^). 

We introduce the operator 

>CBs(tr) = dt + (3.5) 

The following proposition gives the first-order correction to the expression of the 
martingale Mt when 5 is small. 

Proposition 3.1. When 5 is small, we have 

Mt = Qt{Xt) + 0{S^), (3.6) 

where 

Qt{x) = Qi°\x) + 6a(j)t{x'^dlQ[°\x)) +dpQ[^\x), (3.7) 

(x) is deterministic and given by the Black-Scholes formula with constant volatil¬ 
ity a, 

^Bs{^)Qi°\x) = 0, Qt\x) = h{x), (3.8) 

(pt is the random component 


= E 
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(3.9) 




(3.10) 


and (cc) is the deterministic correction 

= a‘^xdx{x^dlQf\x))Dt,T, 

with Dt^T defined by 

r / V \H+- 1 

D,r=ViT-tl = e-(l--) *-4 (3.11) 


The correction solves the problem in (3.16) below. The function 'D{t) derives 
from solving this problem and is: 

^('^) = / ~ u)IC{u)du, 

Jo 

it is discussed in more detail in Lemma A.2 in Appendix A. 

Note that the stochastic volatility process we have introduced in Eq. (3.2) is a 
stationary power-law process. As a consequence of our modeling we have in particular 
that (pt is a Gaussian J^j-measurable process, it reflects the influence of the past on 
the future stochastic volatility path conditioned on the present. We next present 
the proof of Proposition 3.1 and remark that in the analytic framework that we set 
forth, exploiting the “er-martingale decomposition” [27], the cases with H <1/2 and 
H >1/2 can be treated in a uniform way. 

The proof we present below holds for general Gaussian processes Zt with short- 
and long-range correlations, while the expression to the right in Eq. (3.11) is specihc 
to the fOU process. We discuss the general Gaussian case and the general expression 
(B.4) of T>(t) in Appendix B. 

Proof. For any smooth function qt{x), we have by Ito’s formula 

dqt{Xt) = dtqt{Xt)dt+ {xd^qt){Xt)atdW/ + ^{x'^dlqt){Xt)a/dt 
= dlBsi<^t)qtiXt)dt + {xd^qt){Xt)atdW/, 
the last term being a martingale. Therefore, by (3.8), we have 

dQ[°\Xt) = {5aZ^ + y3'(Zf)) {x^dl)Q^°\Xt)dt + dN^°\ (3.12) 

with a martingale, 

divf) = {xd,Qi°^){Xt)atdW/, 

and g^{y) is the function 


9\y) = 2d 


F{Sy) - dy 


FiSy)^ 

(52 


that can bounded uniformly in 6 by 


\g\y)\ < (d||E"||oo + ||A'||L)y". 


Note also that in Eq. (3.12) (and below) we use the notation 
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x=Xt- 





Let (j)t be defined by (3.9). We have 


04 



Z^ds, 


where the martingale ipt is defined by 


V't = 



Z^ds\Ft 


and it is studied in Appendix A. We can write 


Z^{x^dl)Qf\Xt)dt = {x^dl)Q<f\Xt)d^t - {x^dl)Qf\Xt)dct>t. 


(3.13) 


(3.14) 


By Ito’s formula: 

d{<j,t{x^dl)Qf\Xt)) = {x^dl)Q<f\Xt)d^t 

+ {xd4x^dl))Q\°\Xt)atcl)tdW: 

+ ^{x‘^dl{x^dl))Q[°\Xt)(Tt^tdt 
+ [x'^dldt)Qf\Xt)(t)tdt 
+ {xd,{x^dl))Qf\Xt)atd{(l^,W*), 

= {x^dl)Qf\Xt)dcj,t 
+ {xd,{x^dl))Qf\Xt)at<l>tdW: 

+ {5aZ^ + \5^g\Z^)){x'^dl{x^dl))Q‘f\Xt)(i)tdt 

+ {xd,{x^dl))Q^^\Xt)atd{4>,W*),, 

where we have used again £Bs(d)Qt°^(2^) = 0. We have {4>,W*)^ = p{%l),W)^ and 
therefore 

d{<^t{x^dl)Qf\x^)) = -Z^{x^dl)Qf\Xt)dt 

+ {SdZ^ + ^S^g\zi^)){x^dl{x^dl))Q[°\Xt)(l)tdt 

+p{xdx{x'^dl))Q[°^{Xt)atd{i{j, W)^ 

+dNi"\ 

where is a martingale, 

= {xd,{x‘^dl))Q[°\Xt)at4>tdW: + {x^dl)Qf\Xt)d'4^t. 


Therefore: 

d{Qf\x^)+5dcl,t{x^dl)Qf\Xt)) 

= {5‘^d‘^Z^ + ]^5^dg\Z^)){x^dl{x^dl))Qf\Xt)(t)tdt 

+ ^9\Z^){x‘^dl)Qf\Xt)dt + 5dp{xdx{x^dl))Q^°\Xt)<7td{ilj,W)^ 

+dN[°^ + d5dN[^\ (3.15) 
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The deterministic function q [^'^ defined by (3.10) satisfies 


>CBs(o-)Qr^(a:) = -a‘^{xda:{x‘^dlQ[°\x)))0t,T, Qt\^) = 0^ (3-16) 

where 9t,T is such that 


d{'tp,W)^ = 9t,Tdt, 

and it is given by (see Lemma A.l): 

6t,T = [ K.{v — t)dv = f K{v)dv. (3-17) 

Jt Jo 


Applying Ito’s formula 

dQ[^\Xt) = CBs{<Jt)Q[^\Xt)dt + {xd,Q[^^){Xt)atdW: 

= £Bs{d)Q[^\Xt)dt + {5dZ^ + ^g^{zn) ix^dl)Qi^\Xt)dt 
+ {xd,Qi^^){Xt)atdW; 

= -d"^{xda:{x'^dl))Q[°\Xt)d 

+ {SdZi^ + yg^(Zf)) {x^dl)Q\^\Xt)dt + dNi^\ 


f 2) 

where is a martingale, 


dK 


( 2 ) 


{xd,Q[^^)iXt)atdW*. 


Therefore 

d{Qi°\Xt) + Sdcl)t{x^dl)Q[°\Xt) + 6pQi'^\Xt)) = dNt - dRt,T, (3.18) 
where Nt is a martingale, 

Nt= [ dN^°'> + aSdN^ + pSdN^, (3.19) 

and Rt^T is of order 5^: 

Rt,T = £ {d^Z^ + ^dg^Z^)) {x^dl{x^dl))Q^°\Xs)cl,sds 

+ g\Z^) {x^dl) Qf^ {X,)ds + £ dp{xd, {x^dl)) Q(°) {X,)Z^9,,Tds 

+<52 (paZf + ^-pg\Z^)) {x^dl)Q^p{X,)ds. (3.20) 

Then with Qt{x) defined as in Proposition 3.1 we have Qt{x) = h{x) because 
Q^\x) = h{x), (l)T = 0, and Q^\x) = 0. Therefore 

Mt = E[h(XT)\Rt] = E[QT{XT)\Xt] = Qt{Xt)+E[NT - Nt\Tt] +E[i?t,T|^t] 

= Qt{Xt)+E[Rt,T\Xt], (3.21) 

which completes the proof since E[i?t_T|-7t] is of order S^. □ 
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4. Accuracy with European option. In the derivation above we assumed a 
smooth payoff function. Since important classes of payoff functions have non-smooth 
payoff we generalize here the proof to such a class by considering a European option. 
For a European option h{x) = {x — K)+ we have from Eq. 1.41 in [27] 

<«) 

where $ is the cumulative distribution function of the standard normal distribution. 
We can see that h is not smooth so that the hypotheses of Proposition 3.1 are not 
satisfied. However the conclusions of Proposition 3.1 still hold true as we now show. 
Proof. One has to show that Rt^r defined by (3.20) satisfies is of order 

in RP for any p and that the local martingale Nt defined by (3.19) is a martingale 
(up to time T). The problem comes from the fact that the derivatives of q[^\x) blow 
up when t ^ T. However this blow up is not strong as we show below. We first state 
a few properties of the deterministic and random terms that appear in the expression 
of Rt,T- 

- The deterministic function q[^\x) given by (4.1) satisfies 

5,^Q®(x)<C'(l + ^^—^), 

for any 1 < fc < 4, t G [0,T], x G (0,oo), and for some constant C (see Appendix B 
in [25]). 

- The deterministic quantity given by (3.11) satisfies 

Dt,T<C{T-tf+l, 

for any t G [0,T] and for some constant C (see Lemma A.2 and below in Appendix 
A). 

- The deterministic quantity 9t^T defined by (3.17) satisfies 

0t.T<CiT-t)^+^, 

for any t G [0,T] and for some constant C (substitute (2.9-2.10) into (3.17)). 

- The random component (ft defined by (3.9) satisfies 

<Cp{T-t), 

for any t G [0,T] and for some constant Cp for any p > 0 (apply Lemma A.3 in 
Appendix A and use the fact that (ft is Gaussian). 

- The random process satisfies 

E[\Z^\P]i <Cp, 

for any t G [0,r] and for some constant Cp for any p > 0 (use the fact that Zf^ is 
Gaussian, stationary, with mean zero and variance cr^^). 

As a consequence, the deterministic function (x) satisfies 

\d^,Q^t\x)\ < C({T - + (T - t)^+^-4) (1 + x3), 
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for any 1 < A: < 2, f G [0, T], a; G (0, oo), and for some constant C. 

Using (3.20) and the previous estimates we find that, for any p > 0, there exists a 
constant Cp such that 

T 

< Cp5^ / + {T- s)"^ + (T - s)^-^ + (T - sf-ids 

< Cp(52((r - (T - t)"+5), 


for any <5 G (0,1) and t G [0,T], which shows the desired result for Rt,T- 
Moreover, the local martingales in (3.19) are continuous square-integrable mar¬ 
tingales up to time T whose brackets are 

j = 0,l,2, 

AT® = , 

= {{xd^{x^dl))Q^^\Xt)cjt(j)tf 
+ 2p0t.T((xa.(x292))Q(O)(XOu*</,t)((x2a2)Q(O)(XO) 
+{{x^dl)Qf\Xt)felT, 

= (at(xa,QW)(Xt))", 

where the are uniformly bounded with respect to t G [0, T] in for any p, which 
concludes the proof. □ 

5. The implied volatility. We now compute and discuss the implied volatil¬ 
ity associated with the price approximation given in Proposition 3.1. This implied 
volatility is the volatility that when used in the constant volatility Black-Scholes 
pricing formula gives the same price as the approximation, to the order of the approx¬ 
imation. The implied volatility in the context of the European option introduced in 
the previous section is then given by 


It = a + 6 


T-t 


SpDtp 


+ 


log{K/Xt) 


l2{T-t) d{T-ty 


+ 0{S^). 


(5.1) 


The first two terms can be combined and rewritten as (up to terms of order S^): 


u + 5 


T-t 


= E 


r 1 
-T-t 



0{5^). 


(5.2) 


When a(T — t) <C 1 the implied volatility is random and we have (see Lemma 
A.3) and Eq. (A.5) : 


It = (J + SZ. 


r(Er + |) L2 




^ogjK/Xt) - 


(5.3) 


Note that, for H G (0,1/2), the implied volatility blows up at small time-to- 
maturity T — t. Note, moreover that the result above is valid in the asymptotic 
regime ^ <C 1. Indeed, for a being an order one strictly positive quantity the implied 
volatility in Eq. (5.3) is strictly positive for 6 small enough. 

When a{T — t) ^ 1, the quantity Ht.T is of order (T —1)'^+2 and is deterministic 
(by Lemma A.2), while the fluctuations of (f>t are of order (T — t)^ at most and are 
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therefore negligible (by Lemma A.3). As a consequence, when a(T — t) ^ 1, we can 
write the implied volatility as: 



\og{K/Xt) 


(5.4) 


Note that, for H G (1/2,1), the implied volatility blows up at large time-to- 


maturity T — t. We remark that the factors multiplying the square brackets 


in Eqs. (5.3) and (5.4) are slightly modified in the general case when Zt is a general 
Gaussian process, see Appendix B. 

6. A slowf volatility factor. We show in this section that the approach de¬ 
veloped in this paper can be applied to other stochastic volatility models. Here we 
consider the following model 



( 6 . 1 ) 


where E is a smooth, positive-valued function, bounded away from zero, with bounded 


derivatives, and is a rescaled fOU process: 

dZf’^ = S^dWf - SaZl'^dt, 


( 6 . 2 ) 


whose natural time scale is 1/5. It has the form 





(6.3) 


— OO 


Its moving-average integral representation is 



(6.4) 


where K is defined by (2.8). In particular its variance is Uou defined by (2.5), indepen¬ 


dently of 5. This model is therefore characterized by strong but slow fluctuations of 
the volatility. If the price of the risky asset follows the stochastic differential equation 
(3.1), we get a result similar to Proposition 3.1. 

Proposition 6.1. When 5 is small, denoting ao = F{Zq^) and po = F'{Zq^), 
the option price (3.4) is of the form 


Mt = Qt{Xt) + 0{5^^), 


(6.5) 


where 


Qt{x) = Qf\x) + aopQ(t)l[x'^dlQf\x)) + 5^ ppqQ^^\x) , (6.6) 

q|°^(x) is given by the Black-Scholes formula with constant volatility ao, 


t^Bs{cro)Q[°\x) = 0, Qt\x) = h{x), 


(6.7) 


(ff is the random component 



( 6 . 8 ) 


15 








(6.9) 


and (x) is the correction 


= (jlxd^{x'^dlQf\x))Dt,T, 


with Dt^T defined by 


We remark that we indeed can expect the situation with a slow volatility factor to 
behave qualitatively as the situation with small volatility fluctuations in Proposition 
3.1 from the point of view of the effect the medium roughness. This follows since we 
have from self-similarity of fractional Brownian motion that in distribution 



However, we have 



and thus the models (small volatility fluctuations versus slow) differ both in a strong 
sense and in distribution. Moreover, the models have different interpretations from 
the modeling viewpoint with for instance a different skewness mechanism. Note in 
particular that this difference manifests itself in that for fixed Hurst coefficient H the 
magnitude of both the correction and the error terms, deriving in particular from the 
modeling of correlation, have a different scaling in 5 for the two models. For instance 
for small Hurst exponent H we may expect, for given S, the correction (and also the 
error term) to be relatively larger in the case of the slow volatility factor. The random 
correction is of order 6^. More exactly it is a zero-mean Gaussian random variable 
with variance 



+ "]'rfu + 0(<52^+i), (6.11) 


for t G [0,T]. 

Proof. We note that 


(Tt — cTo + — Z^’^) gf, 


where 5 ,^ = - F{Zl'^) - F'{Zl'^){Z^f^ - Z^'^) and therefore 



We have 



which is of order 5^^ ■. 


]E[(Zf’-^ - Zl'^f] = aj^idtf^ + o{S^^). 
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Therefore is bounded in for any p by a quantity of order 5"^^. We can then 
follow the same proof as the one of Proposition 3.1. The term 


dIt= [ {t — u)IC^ {u)du, 

Jo 

is given by 

The variance of the correction (pf is 

=y (y JC^{s — u)ds^ du + J ^y IC^{s — u) — IC^{—u)ds^ du, 


which in turn gives (6.11). □ 

Proceeding as in the case of the small-amplitude stochastic volatiliy model, we 
find that the implied volatility in the context of the European option is given by 


It = cro+po- 


sH 


PPo 


cro 


{T-t) 




\og{K/Xt) 


+ 0 { 5 ^^). ( 6 . 12 ) 


'T-f~ r(if + |)L2^^ 

The first two terms can be combined and rewritten as (up to terms of order S^^)\ 

1 


CTO+PO 


T-t 


= E 


iT-t 


aids\Tt +0{S^^). 


(6.13) 


7. Conclusion. We have presented an analysis of the European option price 
when the volatility is stochastic and has correlations that decay as a fractional power 
of the time offset. The stochastic volatility model is defined in terms of a fractional 
Ornstein Uhlenbeck process with Hurst exponent H and the analysis is carried out 
when the typical amplitude of the volatility fluctuations is relatively small. Two 
situations are differentiated. First the situation when H G (0,1/2) which corresponds 
to a “short-range” dependent process that is rough on short scales with correlations 
that decay very rapidly, faster than linear decay, at the origin. Second the situation 
when H G (1/2,1) so that the correlations decay relatively slowly at large scales 
and then the volatility correlations are not integrable. We use a martingale method 
approach to derive a general expression for the Black-Scholes price covering the two 
cases. In the short-range case the rough behavior on short scales gives rise to an 
implied volatility that diverges as the time to maturity goes to zero. In the long- 
range case the slow decay in the correlations gives a term structure of the implied 
volatility that diverges as time to maturity goes to infinity. The main result we have 
presented is specific in the sense that a particular stochastic volatility model has been 
addressed, however, as we illustrate the framework can be adapted to related models 
as long as some central covariance terms can be computed. We illustrate this by 
considering a model with slow, but order one, volatility fluctuations and derive the 
associated fractional implied volatility term structure. 

Appendix A. Technical lemmas. In this appendix we state and prove a few 
technical lemmas related to some central quantities of interest that are used in the 
derivation of the price in Sections 3 and 5. 
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The martingale ipt is defined for any t G [0,r] by (3.14). It is used in the proof 
of Proposition 3.1 and it has the following properties. 

Lemma A.l. ('i/^t)te[o,T] is a Gaussian square-integrable martingale and 

d ('0, lL)t = ^ y K.{s)ds'^dt, d{'il})^ = (^J K.{s)ds^ dt. (A.l) 

Proof. For t < s, the conditional distribution of given Tt is Gaussian with 
mean 

c* 


E[Z^\Tt] = [ IC{s-u)dWu, 

J —oo 

and deterministic variance given by 

pS — t 

Var(Zf IJ-t) = / lC{ufdu. 

Jo 


(A.2) 


Therefore we have 


f Z^ds+ [ E[Z^\J^t]ds 
Jo Jt 

= f ds f ]C{s — u)dWu f dt f ]C{s — u)dWq 

Jo J — oo Jt J — oo 

= f \ r K{s-u)dt 

J — oo ^ J 0 


’ —OO 

pt , pT 


/C(s — u)dt 


0 J u 


This gives 


p'^ 2 

d = (^ J }C{s — t)ds'^dt, d{'if)^ = (^J }C{s — t)ds'^ dt, 

as stated in the Lemma. □ 

We define the deterministic component 

that appears in Equation (3.10). It has the following properties. 

Lemma A.2. D^ t is a deterministic function ofT—t and it is given by 

Dt,T =T^{T — t), 'D{t) = f {t — u)lC{u)du. 

Jo 

The function V can be written as (3.11) and it has the following behavior: 
For ar <C 1, 


For or » 1, 


V{t) = 


r(Lf + f)a"+i 


)^+5 +o((ar)^+5)). 


H+k 


(A.3) 


(A.4) 


(A.5) 


(A.6) 


Finally, we consider the random process (ft defined by (3.9). 
Lemma A.3. 
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1. 4>t is a zero-mean Gaussian process with variance 


Var((/)t) 



T-t 


IC{s + u)ds^ 


2 

du. 


2. There exists a constant C (that depends on H) such that the variance of (ft 
can be bounded by 


Vax{(ft) <C{T- tf^ A (T - tf. 


3. (ft is approximately equal to {T — t)Zf^ for small T — t: 


E 


(ft 

T-t 



T-t^O 


0 . 


(A.7) 


(A.8) 


Proof. We can express the variance of (ft as: 

VKT{(ft)= [ ds [ ds'Cov{E[Z^\Po],E[Z^\Po]), 
Jo Jo 

which gives the first item since 


Cov(E[zf|J-o],E[zi?|J-o]) 



u)K.{s' 


u)du. 


Furthermore 


Var(0i) < ( J Var(E[Zf 

<(/ IC{u)^duy^^ds^ 


<C{T-tY^ A (T-t) 


which gives the second item of the lemma. 
Similarly, we have 


E 


^ (T^tJo *^®Var(E[zf|J-o] 


and 


/■O 

E[ZY\To] -Z^ = / (/C(s - - /C(-M))dW„, 

J — oo 


so that 


E 




IVT-t 


As s —>■ 0, we have (/C(s-i-v) — /C(v))^dv —>■ 0 by Lebesgue’s dominated convergence 

theorem (remember K € A^), which gives the third item. □ 

Appendix B. Extension to a general stochastic volatility model. In the 

paper, we model the volatility as a bounded function of a fOU process. In fact it is 
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straightforward to extend all the results to a volatility model that is a bounded func¬ 
tion of a stationary Gaussian process whose correlation properties are qualitatively 
similar as the ones of a fOU process. In this appendix we consider the situation when 
the volatility is 

at = a + F{6Zt), (B.l) 

for Zt a stationary Gaussian process with mean zero of the form 

Zt= ! K{t-s)dWs, (B.2) 

J —CO 

where Wt is a standard Brownian motion and K € L^(0,oo) is a general kernel instead 
of the specific kernel (2.8) corresponding to a fOU. Then the Gaussian process Zf has 
mean zero, variance 

}C^{u)du, (B.3) 

and covariance 

pOO 

E[ZtZt+s] = / /C(u)/C(u -I- s)du. 

Jo 

As before (above Proposition 3.1), the function F is assumed to be one-to-one, smooth, 
bounded from below by a constant larger than —ct, with bounded derivatives, and such 
that F(0) = 0 and F'{0) = 1. Proposition 3.1 then holds true, with the function 1) 
defined by 

^('^) = / (''” ~ u)IC{u)du, (B.4) 

JO 

and the implied volatility in the context of the European option is still given by (5.1) 
with Dt^T = T’(T — t). The behavior of the function V is determined by the one of 
the kernel /C and we consider in more detail two cases corresponding respectively to 
long- and short-range correlations: 

1. There exists 0 such that 

/C(t) = czt^~^ (1 -I- 0 ( 1 )) as t ^ 00 . (B.5) 

It H G (1/2,1) this implies that /C is not integrable at infinity and, as we will 
see below (see Lemma B.l), the covariance function of Zt has a tail behavior 
similar to that of a fOU at infinity. In other words, Zt possesses long-range 
correlation properties, and the implied volatility has the same form (5.4) as 
in the case of a fOU with Hurst index H, with czT(H — l/2)/r(iL -|- 3/2) 
instead of l/[ar(iL -1-3/2)]. 

2. There exists dz ^ 0 such that 

/C(t) = dzt^~^ (1 -I- 0 ( 1 )) as t —>■ 0. (B.6) 

It FI G (0,1/2) this implies that K. is singular at zero and, as we will see 
below (see Lemma B.2), the covariance function of Zt has a behavior similar 
to that of a fOU at zero. In other words, Zt possesses short-range correlation 
properties, and the implied volatility has the same form (5.3) as in the case of 
a fOU with Hurst index H, however, with dztt{F[-\-l/2)/T{FIFt)/2) replacing 
l/r(iL-f 5/2). 
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Lemma B.l. We assume (B.5). 

1. If H G (1/2,1), then the covariance function of Zt satisfies 

E[ZtZt+s] = + o(l)) as oo, 


with 


_ , r(2-2g)r(g-i) 2 

^ r(|-iL) ^^2sin(7riJ)r(2iJ-!)■ 

2. If H € (1/2,1), then the function V^t) defined by (B.4) satisfies 

(1 + 0 ( 1 )) asT^oo. 

1 (H + 2 ) 


(B.7) 


(B.8) 


(B.9) 


If Zt is the fOU process (2.4), we have cz = l/[ar(iJ— 1/2)]. In this case, we can 
check that kz = a~^ /[2sui{TrH)T(2H — 1)] = /T{2H — 1), which confirms 

that (B.7-B.8) give (2.12), while (B.9) gives (A.6). 

Proof We denote 

poo poo 

C{s) = E[ZtZt+s\ = / /C(m)/C('« + s)^^ and C(s) = {u + s)^~i du. 

Jo Jo 


We can check that C(s) = kzs^^ ^ with 


kz = 


f 


.H-i 


{1 + uW ^du = c] 


r (2 - 2i7)r(i7 - i) 


r(| 


-H) 


We now show that C{s) —C(s) goes to zero as s —00 faster than s^^ Let £ G ( 0 , 1 ). 
There exists such that |/C(f)t“^+^/^ — cz\ < £ for any t > . We have for any 

s > S^-. 

/•S' 

s^“^^|C(s) — C(s)| < / |/C(m)/C(m + s) — C 2 U^“ 5 (u- p s)^-| 

JO 

poo 

/ |/C(u)||/C(m + s) — C2(u + s)^“^|(ilt 

JS^ 

pOO 

_^g2-2H / \cz\{u + s)^~^IC{u) - CzU^~^du 
JS'^ 

^ ^2-2H / ^}C{u)}C{u + s) — C%U^~i {u + s)^~^\du 

Jo 

poo 

_l_g 2 - 2 ff£( 2 |^^| _l_ £) / (^u + s)^~iu^~idu. 

Jo 

As s —>■ 00 the first term of the right-hand side goes to zero by Lebesgue dominated 
convergence theorem because (2 — 2H) + {H — 3/2) < 0 . This gives 


limsups^ ^'^|C(s) — C(s) 

s—^■oo 



idu. 


Since this holds true for any £ G (0,1), this proves (B.7). 
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We denote 


Vir) 



. JT _ ^ 

ijczu ^du, 


which is given by 


V{t) 


Cz H+i- 

_T ^ 2 


Cz 


T{H - 
T{H + 



Let e G (0,1). As mentioned above, there exists 5"^ such that \K.{t)t < e 

for any t > S^. We have then for any r > S'®: 

/■S' 

— 'D{t)\ < T~^~^ / {t — u)\}C{u) — czu^~^du 
Jo 

+T~^~2£ f (T — u)u^~idu. 

Js<^ 

As T —>■ oo the first term of the right-hand side goes to zero by Lebesgue dominated 
convergence theorem because 1/2) -|- 1 < 0. This gives 


limsupr ^ —'b{T)\ < e f {1 — u)u^ ^du. 

T—¥00 JO 


Since this holds true for any e G (0,1), this proves (B.9). □ 

Lemma B.2. We assume (B.6). 

1. If H G (0,1/2) and if 1C satisfies the two technical conditions: 

(CB.2.1) fC is integrable and Lipschitz on (l,oo). 

(CB.2.2) There exist functions ki{t) and k 2 {s) such that for all t,s G (0,1) we 
have \IC{t -b s) — IC{t)\ < ki{t)k 2 {s), where iC{t) = lC{t) — dzt^~^^^, 
ki G -L^(0,1), and liins^o s~^k 2 {s) = 0. 

Then the covariance function of Zt satisfies 

^[ZtZt+s] = cr\ — qzs"^^ oi^s^^) as s —>■ 0, (B.IO) 


with 


4 r(g + i)^ 

2 r( 2 iL -b 1) sin(7riJ) ’ 



2. For any H G (0,1), the function 'D(t) defined by (B.4) satisfies 


V{t) 


+ asr^O. 


(B.ll) 

(B.12) 


(B.13) 


The condition (CB.2.1) gives some control of /C away from the origin, and this 
specific condition can be relaxed. The necessary condition is that 


„-2ff 


P°° 9 

J (/C(u-b s) —/C(m)) du 
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goes to zero as s —>■ 0 (see the proof below). 

The condition (CB.2.2) means that the remainder lC{t) should be small enough 
near the origin. A sufficient condition for (CB.2.2) is that 1C is a-Holder continuous 
over (0, 2) for some a > H. Then (CB.2.2) is fulfilled with ki{t) = c and fc 2 (s) = kaS°‘, 
for some constant c. 

If Zt is the fOU process (2.4), then 1C is integrable and Lipschitz on (1, oo) and we 
have dz = l/r(i7 + 1/2) and iC{t) = —[a/r{H + 1/2)] f^(t — which is 

{H + l/2)-H61der continuous over (0, 2): \iC{t + s) — ^(t)| < [2a/r(iJ + 3/2)]s^+^/^. 
In this case we can check that qz = l/[2r(2iJ + l)sin( 7 riJ)] = /T{2H + 1), 

moreover we have then cr| = o-^^, which confirms that (B.10-B.12) give (2.11), while 
(B.13) gives (A.5). 

Proof. We can write 

E[ZtZt+s] = 4 - Q{s), Q{s) = ^E[{Zt+s - Ztf]. 

We have Q(s) = Qi(s) + 62 ( 5 ) with 

Qi(s) = - / (/C(u-l-s) -/C(u))^du, Q 2 (s) = - / /C(u)^du. 

^ Jo ^ Jo 

The idea is to approximate these two functions by their versions when dzt^~^^^ 
replaces /C(t). We denote 


poo 

Qi(s) = 4 / ((u + s)^-i - u^-i^du, 

^ Jo 

We can check that Qi(s) + Q 2 (s) = qzs^^ with 
qz = -^j^ {u^ ^ -{1+uf du+^ 


62(5) = ^ ^ du . 

^ Jo 


u'^H-ld ^ 4 + 2 )^ 

2 r(2iJ + 1) Sin(7ri/) ■ 


We now show that Qi(s) — Qi(s) goes to zero as s —>■ 0 faster than We have 

2s-""|Qi(s)-Qi(s)| 

f {lC{u + s) — lC{u)Y — d\i^{u + s)^~i — u^~i)^du 
Jo 

poo 

[{u + s)^~i — u^~iYdu + 


< s 


-2H 


+ 5 - 2^4 


< 2\dz\s / ((u + s) 


\H-h . R-1\2^ 


ll/2r 


(/C(u + s) —/C(u)) du 
1/2 


) du / [K(u + s) — ^(w)) du 


+s I [}C{u + s) — }C{u)ydu + s ^^4 / ((“ + 'S)’^ ^ — u^ du 


/ OO 

|/C(u + s) — /C(rt)| {\lC{u + s)| + \IC{u)\)du 


<2\dz\ I {{u+lf-^ -u^-y^du 


1/2 r 


,-2i? 


[lC{u + s) — IC{u)ydu 


-,1/2 


+s 


-2H 


(^^{u + s) — ^{u)Ydu + d^z I (('“ + 1)’^ — u^ du 

Jljs 


/'OO 

+2s^~^^Ltc / \IC{u)\du, 
Jo 
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where Ljc is the Lipschitz constant of /C over (1, oo). As s —^ 0 the third term of the 
right-hand side goes to zero because the integral is convergent and the fourth term 
goes to zero because 1 — 2H > 0. The first and second terms go to zero because 

/ {lC{u + s) — IC{u))'^du < k2{s)^ / ki{u)^du, 

Jo Jo 

ki € L^(0,1), and s~^k 2 {s) —>■ 0 as s —0. Therefore 

lim - Qi(s)| = 0. 

s^O ' ' 


We now show that Q 2 {s) — Q 2 {s) goes to zero as s —0 faster than . Let e € (0,1). 
There exists 5"^ such that — dz\ < e for any t < S'^. We have for any 

s < S‘^: 


2s-"" Q2 (s) - Q2(s) <s- 


-2H 


K.{u) — dnu^ ^\{\K.{u)\+ \dz\u^ ^)du 


<s e{2\dz\ + e) f ^du< 

Jo 


(2|dz| -I- s)e 
Tff ■ 


Since this holds true for any e G (0,1), we have 


lims 2"|Q2 (s) - Q2(s)| = 0 , 


which completes the proof of (B.IO). 
We denote 


=dz (r - u)u 
Jo 


H- 


2 du. 


which is given by 

V{t) = 


(iL + i)(Lr + |) 




= dz 


r(7L + |) 


2)H+i 


Let £ G (0,1). There exists such that \K.{t)t ^+1/2 — dz} < e for any t < S^. We 
have for any t < S^: 


-H-4 


= T>(t) -V{t) \ < 




(r — u)'. 


.H-i 


^du. 


This gives 


lim sup r — 'D{t)\ < e / (1 — ^du. 

r-s-O Jo 


Since this holds true for any e G (0,1), this proves (B.13). □ 
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